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Molecular dynamics simulations are performed to study the temperature-dependent dynamics and 
structures of the hydration shells of elastin-like and collagen-like peptides. For both model peptides, 
it is consistently observed that, upon cooling, the mechanisms for water dynamics continuously 
change from small-step diffusive motion to large-step jump motion, the temperature dependence 
of water dynamics shows a weak crossover from fragile behavior to strong behavior, and the order 
of the hydrogen-bond network increases. The temperature of the weak crossover from fragile to 
strong behavior is found to coincide with the temperature at which maximum possible order of 
the hydrogen-bond network is reached so that the structure becomes temperature independent. 
In the strong regime, the temperature dependence of water translation and rotational dynamics 
is characterized by an activation energy of Ea ~ 0.43 eV, consistent with results from previous 
dielectric spectroscopy (DS) and nuclear magnetic resonance (NMR) studies on protein hydration 
waters. At these temperatures, a distorted 7r-flip motion about the twofold molecular symmetry 
axes, i.e., a water-specific P process is an important aspect of water dynamics, at least at the water- 
peptide interfaces. In addition, it is shown that the hydration waters exhibit pronounced dynamical 
heterogeneities, which can be traced back to a strong slowdown of water motion in the immediate 
vicinity of peptide molecules due to formation of water-peptide hydrogen bonds. 



I. INTRODUCTION 

A number of anomalies distinguish water from other 
liquids. ^'^ For example, it was argued that supercooled 
water exhibits a fragile-to-strong transition (EST) in the 
vicinity of = 228 K^ii Here, FST means that super- 
cooled water behaves as a fragile liquid above T^, i.e., 
the temperature dependence of viscosity deviates from 
an Arrhenius law, while water is a strong liquid below Tx 
and, hence, an Arrhenius law is obeyed.^ It was put for- 
ward that the proposed FST is related to a hypothesized 
liquid-liquid critical point at elevated pressure, which ter- 
minates a phase transition line between a low-density liq- 
uid (LDL) and a high-density liquid (HDL) at low and 
high temperatures, respectively.^ However, direct obser- 
vation of the proposed FST is not possible for bulk water 
because of inevitable crystallization?^ 

By contrast, crystallization can be suppressed in con- 
finement so that the dynamics of supercooled waters are 
accessible down to the glass transition. Confined and 
interfacial waters are of enormous importance for biolog- 
ical, geological, and technological processes. While some 
experimental studies on such waters observed crossovers 
in the temperature dependence of a correlation time r 
at T « 224 K"^'* or T « 200 others provided evi- 

dence for an absence of any crossover jii Several work- 
ers took the crossover at ca. 224 K as indication for a 
FST and related it to a liquid-liquid transition of con- 
fined waters F^iSiii^ Other workers challenged this conclu- 
sion and argued that the Arrhenius processes found below 
the crossovers are secondary (/3) relaxationS )^°i^^ such as 
the Johari-Goldstein (JG) (3 process)^ rather than the 
structural {a) relaxation, which is difficult to observe at 
these temperatures 

Protein hydration waters are prominent examples of 



interfacial waters. Due to an interplay of protein and wa- 
ter dynamics, existence of a hydration shell is essential 
for the biological functions of proteins.-^ Despite appre- 
ciable progress in recent years, the nature of this cou- 
pling is still controversially discussed. It was argued 
that the protein dynamics are "slaved" by the water 
dynamicsj ^^'^^ However, it is not clear which relaxation 
processes of protein (local or global) and water (a or 
13) are related i-^°i^^i^^'^°i^^i^^ Valuable insights have been 
gained by characterizing the temperature dependence of 
protein and water dynamics, in particular during the so- 
called dynamical transition of proteins at ca. 200 K )^°i^^ 
which was reported to be associated with a freezing of bi- 
ological functions.^'' Interestingly, it was postulated that 
the dynamical transition of proteins is related to a FST 
of the hydration waters4i2^ 

A fundamental understanding of the coupling of pro- 
tein and water dynamics requires, first, assignment of the 
various relaxation processes to protein and water, respec- 
tively, and, second, identification of the processes as a or 
(3 relaxations. Recently, we demonstrated that ^H NMR 
is a powerful tool for this purposei^ For elastin and col- 
lagen, two important proteins of the connective tissue, 
we found that, when the temperature is decreased, the 
mechanism for the reorientation of the hydration waters 
changes from an isotropic motion to an anisotropic mo- 
tion, which is dominated by large-angle rotational jumps. 
The mechanism for the low-temperature motion is con- 
sistent with neither the a process nor the JG (3 process 
of molecular glass-forming liquids^^iS^ implying the ex- 
istence of a water-specific (3 process. 

Providing access to the complete microscopic infor- 
mation, molecular dynamics (MD) simulations enable 
straightforward assignment and identification of the var- 
ious relaxation processes. A huge amount of simulation 
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works focused on the behaviors of hydrated proteins un- 
der physiological conditions^^i e.g., the mechanisms for 
hydrogen-bond (HB) breaking were shown to differ for 
bulk and hydration waters?^ Computational studies on 
the temperature-dependent properties of hydrated pro- 
teins are more rare .^^i^" '^'^'^'^'^'^ It was found that 
the protein dynamical transition is triggered by the re- 
laxation of the protein-water HB network, which in turn 
is driven by water translational diffusion i'^^i^^i'^^i^'* More- 
over, the protein dynamic transition was found to occur 
at the temperature of both a dynamic crossover in the 
diffusivity of hydration water and a maximum of specific 
heat, suggesting a relation to the hypothesized liquid- 
liquid critical point of water and to an evolution of water 
from HDL to LDL.— 122 These results imply that water 
plays an important role for protein dynamics. 

We use MD simulations to investigate the temperature- 
dependent interplay of protein and water dynamics for 
hydrated elastin and hydrated collagen. In doing so, 
we exploit that peptides (VPGVG)„ have proven suc- 
cessful models of ideal elastin^^i^Si^ while triple-helical 
molecules [3(P0G)„]'^i model ideal collagen^^sdi The 
choice of the proteins elastin and collagen is motivated 
by both their biological importance and their use in our 
previous ^H NMR study^ so that we expect synergetic 
effects from a combined experimental and theoretical ap- 
proach. In the present contribution, we focus on the 
temperature-dependent translational and rotational mo- 
tion of the hydration waters to improve our knowledge 
about the dynamics of supercooled bulk water and the 
role of water for the biological functions of proteins. The 
protein dynamics and the interplay of the components 
will be investigated in future work. 



II. METHODS 

We study three models with the same hydration level 
h = 0.3 (g water/1 g peptide). The smaller model of hy- 
drated elastin (El) is comprised of 1 peptide (VPGVG)5o 
and 342 H2O, while the larger model (E8) contains 8 pep- 
tides (VPGVG)5o and 2732 H2O. The model of hydrated 
collagen (C5) is comprised of 5 triple-helical molecules 
[3(POG)io] and 705 H2O. 

MD simulations were performed using the GROMACS 
software package4^ The GROMOS96 43a2 force field 
was utilized^^ Moreover, the SPG model of water was 
employed;^ as recommended for GROMOS force fields. 
We applied periodic boundary conditions and a time step 
of 2fs. The nonbonded interactions were calculated uti- 
lizing a cutoff distance of 1.2 nm. The PME technique^ 
was used to treat the Coulombic interactions and the 
LINCS^** and SETTLE*^ algorithms were applied to con- 
strain the bonds of peptide and water, respectively. Prior 
to data acquisition, the systems were equilibrated in sim- 
ulations at constant N, P, and T, using the Rahman- 
Parrinello barostat^'^ and the Nose-Hoover thermostat. 
These equilibration runs, which spanned up to 100 ns for 
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FIG. 1; Pair distribution functions (a) goo and (b) goH for 
E8 at the indicated temperatures. 



El and up to 40 ns for E8, allowed us to adjust the den- 
sities p{T), which increase linearly by 5-7% upon cooHng 
in the studied temperature ranges. The sizes of the cubic 
simulation boxes amounted to about 3.3 nm and 6.7 nm 
for El and E8, respectively. In the case of G5, the atoms 
of the triple helices [3(POG)io] were restrained to their 
positions in the crystal structure, suppressing diffusive 
motion of the peptides. To model a collagen microfib- 
ril, five parallel triple-helical molecules were pentagonally 
grouped, where the distance between the centers of adja- 
cent molecules was ca. 1 nm. The temperature-dependent 
spatial distribution of the water molecules about the 
triple helices was adjusted in equilibration runs at con- 
stant N, P, and T, which spanned up to 40 ns. For C5, 
the sizes of the simulation box were « 8.4 nm and 
Ly Lz ~ 3.4 nm. For all systems, the subsequent pro- 
duction runs were performed in the canonical ensemble, 
i.e., at constant N, V, and T, employing the Nose- Hoover 
thermostat. The trajectories were saved every 1 ps for 
later analysis. 



III. RESULTS 
A. Structure of hydration waters 

To study the structure of the hydration shells, we show 
the intermolecular pair distribution functions goo ("r) and 
9oh{t) for the water oxygens and water hydrogens of E8 
in Fig. [TJ It is evident that the features of goo and 
goH become sharper and an oscillatory behavior devel- 
ops when the temperature is decreased. The positions of 
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FIG. 2: Temperature-dependent probability p of finding a 
water molecule that participates in exactly 4 HB for El and 
E8. 



local maxima and minima, which are essentially indepen- 
dent of temperature, are in good agreement with those in 
previous experimental and computational studies on bulk 
water, including results for SPC bulk water. Specifi- 
cally, the first maximum of goo and goH is located at 
roo — 2.8 A and roH = 1.8 A, respectively. The first max- 
imum of goH results from oxygen-hydrogen pairs forming 
a HB. These findings show that a defined HB network 
develops at the water-peptide interface upon cooling. In 
addition, we observe enhanced values gooij") > 1 and 
goHij) > 1 up to r ~ 25 A, indicating that repulsion of 
water from the hydrophobic cores of the coiled elastin- 
like peptides leads to an agglomeration of water between 
the coils. For none of the models, we find evidence for 
an onset of crystallization at the studied temperatures. 



The temperature-dependent structure of the HB net- 
works can be investigated in more detail, when we de- 
termine the probability p that a water molecule partici- 
pates in exactly 4 HB. Exploiting that gonif) exhibits a 
pronounced first minimum, we define that pairs of oxy- 
gen and hydrogen atoms form a HB when their inter- 
atomic distance is smaller than 2.4 A. In doing so, we 
include water-water and water-peptide HB. In Fig. [21 we 
see that p continuously increases upon cooling until it 
levels off at about 190 K, at least for El. These findings 
suggest that the HB network adopts energetically more 
favorable configurations upon cooling and, eventually, it 
finds one of the best possible configurations under the 
given constraints provided by the water-peptide inter- 
face. The absence of a discontinuity implies that a first 
order phase transition, e.g., between LDL and HDL wa- 
ters, does not occur at ambient pressure. Due to limited 
computer power, it was not possible to equilibrate the 
structure of E8 below 200 K so as to ascertain whether p 
saturates for this model as well. 
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FIG. 3: (a) Incoherent intermediate scattering functions 
So{(l,t) (g = 27r/2.8A) for the oxygen atoms of water in C5 
and El. (b) Orientational correlation functions for the O-H 
bonds of water in C5 and El. In both panels, the results for 
C5 and El are shown as symbols and lines, respectively. For 
C5, the temperatures are indicated. For El, the temperatures 
are 310, 260, 230, 200, 180, and 160 K (from left to right) 



B. Dynamics of hydration waters 

To investigate the temperature-dependent dynamical 
behaviors of the hydration waters, we exploit that the 
incoherent intermediate scattering function 



Siq,t) = (cos{q. [r(to + i) -r(to)]}) 



(1) 



provides us with information about translational motion, 
while the orientational autocorrelation function 



F2(i) = i(3[e(to + t)-e(io) 



-1) 



(2) 



yields insights into rotational motion. Specifically, the 
scattering function S depends on the atomic translational 
displacements [r{to +t) —r{to)] during the time interval 
t, where the absolute value of the scattering vector (/= |q| 
determines the length scale on which dynamics is probed. 
We calculate the scattering functions of the water oxygen 
atoms. So, and the water hydrogen atoms, Sh, for val- 
ues of q corresponding to the respective intermolecular 
interatomic distances. However, we determined that our 
conclusions do not depend on these specific choices of q. 
The orientational correlation function F2 depends on the 
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angular displacements during the time interval t, more 
precisely, on the value of \e{to + t) ■e{to)\- Here, the unit 
vector e(to) describes the orientation of an 0-H bond at 
a time to. In recent NMR work,^^ F2 was determined 
experimentally for the 0-D bonds of heavy water in the 
hydration shells of elastin and collagen. Throughout this 
paper, the brackets (. . .) denote the average over various 
time origins to and over all atoms or bonds belonging to 
the considered atomic or bond species. 

In Fig. [21 we compare temperature-dependent data Sq 
and F2 for El and C5. It is evident that both these func- 
tions exhibit strongly nonexpontial decays, which shift 
to longer times upon cooling. Despite some quantitative 
differences, So and F2 decrease on similar time scales and 
the behaviors of water near elastin-like and collagen-like 
peptides are comparable. Consistently, the water dynam- 
ics of hydrated elastin and collagen samples were found 
to be similar in our previous NMR study,^- For the 
higher of the studied temperatures, the observation of 
complete decays indicates that the a relaxation of hy- 
dration water is probed. Specifically, the results for So 
show that all water molecules move at least one inter- 
molecular distance on a time scale of a few nanoseconds 
at 260 K. Nonexponential a relaxation is a characteris- 
tic feature of supercooled liquids approaching their glass 
transitions. While a stretched exponential function well 
describes the time dependence of the a relaxation in these 
materials, such function does not yield a reasonable in- 
terpolation in our case. We will return to this point in 
Sec lmPl 

To study the slowdown of water dynamics upon cool- 
ing, we determine the translational correlation times tq 
and th and the rotational correlation times th from 
the corresponding incoherent scattering and orientational 
correlation functions. Figure [3] shows the temperature 
dependence of these time constants for El and C5. We 
see that th follows an Arrhenius law with comparable 
activation energies Ea = 0.40 — 0.43 eV for both models. 
Such activation energies, which are typical of breaking 
of 2 HB, were reported for water dynamics in various 
types of confinement at sufhciently low temperatures.^*^ 
The temperature dependence is weaker for tq and at 
high temperatures. Thus, there are deviations from an 
Arrhenius law, which are typical of fragile supercooled 
liquids like weakly supercooled bulk water 1^ Here, the dif- 
ferent temperature dependences of To and th imply that 
diverse mechanisms for water dynamics exist, including 
anisotropic motions. While tq reflects the translational 
motion of the water molecule, the value of t// depends on 
both translational and rotational motions. For example, 
TT flips about the twofold molecular symmetry axis would 
affect Th, but not tq, see Sec. IIII CI 

Next, we discuss the atomic mean-square displace- 
ments (MSD) 




(3) 



Figure shows the temperature-dependent MSD of the 
water oxygens in El, which is a measure for the center-of- 
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FIG. 4: Correlation times of water dynamics in (a) El and 
(b) C5. In both panels, we compare translational correla- 
tion times TO and th for the oxygen atoms from So{<l,t) 
{q — 27r/2.8A) and for the hydrogen atoms from SH{q,t) 
(q — 271/2.4: A), respectively, with rotational correlation times 
tr from F2{t). All shown values are 1/e-decay times, e.g., 
F2{t) — e~^. The dashed lines are Arrhenius fits to th- 



mass motion of water. In addition to a ballistic regime, 
which is expected for short times outside the studied 
time window, three regimes exist for the translational 
motion of water. At short times and low temperatures, 
we observe a plateau regime, indicating that the water 
molecules are trapped in cages formed by neighboring 
water and peptide molecules. This cage effect is typical 
of supercooled liquids approaching their glass transitions. 
At longer times, the water molecules start to escape from 
these cages and, hence, the MSD increases. Subsequent 
to the plateau, we find a regime of sublinear diffusion, 
cxi*^-^^, reflecting the effect that the presence of large 
peptide chains prevents water from free diffusion. In har- 
mony with these results, a recent study on the hydration 
water of RNase A found that a plateau of the MSD devel- 
ops between a regime of ballistic motion and a regime of 
sublinear diffusion when the temperature is decreased.— 
Finally, at high temperatures and long times, there is a 
crossover to linear diffusion, oc t, indicating that the 
spatial confinement due to the peptide becomes unim- 
portant at sufficiently large length scales. These three 
regimes are also observed for water diffusion in C5 and E8 
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FIG. 5: MSD for the oxygen atoms of water in El at various 
temperatures. The hues are power laws {r'^{t)) oc and 
(r-^(t)) oc t"'''^, respectively. 



and, hence, their existence is not an artifact of the small 
size of model El. For El, a confinement effect results 
since the water molecules do not only interact with the 
single peptide chain in the simulation box, but, because 
of the use of periodic boundary conditions, also with all 
its periodic images. The models C5 and E8 contain sev- 
eral peptide molecules so that there are also confinement 
regions within the simulation box. 

When the temperature is decreased, the plateau regime 
extends, indicating that the water molecules become 
trapped for longer and longer times in their local cages. 
Despite this strong delay of water transport, we do not 
find evidence that the diffusion of the hydration wa- 
ter ceases, but it rather exits the time window of the 
simulation in a continuous manner. The relations be- 
tween water diffusion and protein dynamics will be dis- 
cussed in future work. While a self-diffusion coefficient of 
D — 5.6 ■ 10^^" m^/s results from the linear regime of the 
MSD for the hydration water of El at 300 K, a value of 
_D«4.2 • 10~^m^/s was reported for SPC bulk water at 
298Ki^i^ Thus, the presence of the elastin-like peptide 
leads to a slowdown of water diffusion by a factor of 7—8 
at ambient temperatures. 



C. Mechanisms for the dynamics of hydration 
waters 

The mechanisms for the translational motion of water 
can be investigated on the basis of the van Hove self- 
correlation function of the water oxygens 

Go{r,t) = {6[r-\r{to+t)-r{to)\]) (4) 

It measures the probability that a water oxygen moves a 
distance r in a time interval t. In simple cases, the van 
Hove correlation function is a Gaussian. Figure [6] shows 
Go{r] t) of El for various temperatures and comparable 
time intervals, which correspond to late stages of the a 
relaxation, see Fig. [31 It is evident that the Gaussian 
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FIG. 6: Self part of the van Hove correlation function Go{r; t) 
for the oxygen atoms of water in El. For all studied tempera- 
tures, we use comparable time intervals t « lOr^. The dashed 
line marks the interatomic oxygen-oxygen distance roo- 
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FIG. 7: Probability distributions poH{'y',t) of finding an an- 
gular displacement 7 during a time interval t for the O-H 
bonds of water in El. For all studied temperatures, we use 
comparable time intervals t ^tr. The dashed line marks the 
tetrahedral angle. 



approximation is not obeyed at any of the studied tem- 
peratures. Upon cooling, Go{r;t) develops a multi-peak 
structure, where the positions of the secondary maxima 
coincide with the positions of the peaks of gooi^), see 
Fig- HI In particular, a well resolved secondary maximum 
exists at r^rpo- These observations show that, when 
the temperature is decreased, a HE network with defined 
sites develops and the energy barriers between these sites 
start to govern water dynamics, leading to a jump mo- 
tion of the water molecules. Here, the term jump means 
that long periods of vibrational motion about the sites 
are interrupted by short periods necessary to cross the 
saddles that separate the sites in the HE network. 

Valuable information about the mechanisms for the ro- 
tational motion of water can be obtained when we deter- 
mine the probability distribution ^5(7; t) of finding an an- 
gular displacement 7 in a time interval t. In other words, 
p(7; t) describes the probability that the orientations of 
a given bond at the beginning and at the end of the time 
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FIG. 8: Comparison of the probability distributions ps{'y;t) 
and Pp{'y;t) providing information about the temperature- 
dependent anisotropy of water reorientation in El. While ps 
describes the angular displacements of the two-fold symme- 
try axes of the water molecules, pp characterizes the angular 
displacements of a vector perpendicular to the plane defined 
by the three atoms of the respective molecule. For both prob- 
ability distributions, we use comparable time intervals t «r_R 
at the indicated temperatures. 



interval differ by an angle 7. First, we study poHi"f',t) 
for the 0-H bonds of water. In Fig. [3 we display poH 
of El for various temperatures and t m. Upon cool- 
ing, there is a continuous crossover from a one-peak to 
a two-peak signature, indicating that the rotational dy- 
namics of water evolves from small-angle to large-angle 
motion. The position of the secondary maximum shows 
that the O-H bonds tend to jump about angles close to 
the tetrahedral angle at low temperatures. 

To obtain further insights into the nature of the ro- 
tational motion, we investigate the reorientation of two 
orthogonal vectors that are linked to the molecular frame 
of water. Specifically, we calculate ps{'^;t) and pp{'~f;t), 
which describe the rotational motions of the two-fold 
symmetry axis of the water molecule and of the normal 
vector of the plane formed by the atoms of the water 
molecule, respectively. In Fig. [H we show both quan- 
tities for El at various temperatures using t k, t^. We 
see that a secondary maximum exists for pp, but not 
for Ps- The latter finding means that the reorientation 
of the two-fold axis is small during the jump motion of 
most molecules, i.e., this axis essentially coincides with 
the axis of rotation. The position of the secondary max- 
imum of Pp indicates that the normal vector performs 




FIG. 9: Separate analyses of water dynamics in E8 for water 
molecules forming different numbers rip of hydrogen bonds 
with proteins: (a) Incoherent intermediate scattering func- 
tions Soi<l,t) (g = 27r/2.8 A) for the oxygen atoms of water at 
230 K and (b) probability distribution Pp{'^,t) for t = 10ns at 
200 K. 



jumps by about 180°. Together, these results show that 
distorted tt flips about the two-fold symmetry axes are 
an important aspect of water dynamics at sufficiently low 
temperatures, at least in our models. 



D. Effects resulting from the peptide interfaces 

Finally, we study whether water dynamics changes in 
the vicinity of elastin-like peptides. For this purpose, 
we distinguish between water molecules forming differ- 
ent numbers of HB to the peptides in our analyses. In 
Fig. [21 we present the outcome of such separate analyses 
for the example of E8. In panel (a), we compare the inco- 
herent scattering functions Soiq,t) of the different water 
species at 230 K. We see that the decay resulting from 
water molecules that form 2 HB with peptides is about 
three orders of magnitude slower than that resulting from 
water molecules forming no HB with peptides. Thus, as 
a consequence of HB formation, the immediate vicinity of 
a peptide leads to a strong slowdown of the water dynam- 
ics. These pronounced dynamical heterogeneities result 
in a pronounced non-exponentiality of the translational 
and rotational correlation functions for the whole ensem- 
ble of water molecules, see Fig. |31 In particular, they re- 
sult in deviations from a stretched-exponential behavior. 
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which is usually found for supercooled liquids. Specifi- 
cally, a long-time tail results from water molecules that 
are more tightly bound to a peptide through two HB. 
Comparable results are obtained for El and C5. 

In panel (b), the heterogeneity of water dynamics is 
demonstrated for the rotational motion on the basis of 
Pp{'y;t) at 200 K. For water molecules with Up — 0, 
we see that the distribution is close to the final state 
p{'^,t —> oo) oc sin7, indicating that the normal vectors 
have been isotropically redistributed on the unit sphere 
during the used time interval t = 10 ns. By contrast, for 
water molecules with rip = 2, the intensity is cumulated 
in two peaks at small and large angles, respectively, in- 
dicating that most molecules in the immediate vicinity 
of a peptide have only performed distorted tt flips during 
this time interval. 



IV. DISCUSSION AND CONCLUSION 

We have investigated the temperature-dependent wa- 
ter dynamics in the hydration shells of elastin-like and 
collagen-like peptides using MD simulations. Analysis 
of the simulation data shows that the water dynamics 
hardly depend on the type of the peptide. For all stud- 
ied models, water dynamics is strongly slowed down in 
the immediate vicinity of the peptides, in particular at 
low temperatures, because of formation of HB between 
water and peptide molecules, leading to pronounced non- 
exponential correlation functions of water motion. The 
mechanisms for both translational and rotational water 
dynamics change from small-step (diffusive) motion to 
large-step (jump) motion upon cooling. Concerning the 
rotational motion, we have found that distorted tt flips of 
the water molecules about their twofold symmetry axes 
become important when the temperature is decreased, at 
least for molecules at the water-peptide interfaces, i.e., a 
water-speciflc /3 process evolves. While this /3 process 
of water is associated with large-angle jumps, the JG 
/3 process results from small-angle motion^ suggesting 
that both secondary relaxations are not directly related. 
In harmony with the present results, the hydration wa- 
ters of lysozyme were found to perform large-angle jumps 
below ambient temperatures 

Furthermore, we have observed complex water self dif- 
fusion, which is characterized by four time regimes at suf- 
flciently low temperatures: a ballistic regime, a plateau 
regime, a regime of sublinear diffusion, and a regime of 
linear diffusion. When the temperature is decreased, 
the plateau regime extends, indicating that the water 
molecules become trapped for longer and longer times in 
cages formed by their neighbors. A regime of sublinear 
diffusion regime was also observed for the hydration wa- 
ters of RNase A4i It results because the larger peptide 
chains form some conflnement for the water molecules 
and, thus, hamper free diffusion, which is only observed 
at sufficiently large time and length scales. The existence 
of pronounced sublinear diffusion means that determi- 



nation of reliable diffusion coefficients from MD simula- 
tions requires relatively long equilibration and produc- 
tion runs. 

For the temperature dependence of various correlation 
times, e.g., of the rotational correlation time, we have 
found a weak crossover from a non-Arrhenius behavior 
at high temperatures to an Arrhenius behavior at low 
temperatures, which occurs at t« 1 ns. A weak crossover 
was also reported for the water self-diffusion coefficient in 
simulation work using the TIP5P model of water together 
with the GROMOS model of lysozyme. However, the 
present findings are contrary to computational results 
for TIP4P-Ew water in the hydration sheU of OPLS-AA 
lysozyme)^ Despite comparable correlation times of a 
few nanoseconds at the crossover in the present and pre- 
vious studies, a sharp kink toward an Arrhenius law char- 
acterized by a low activation energy of £^Q«0.15eV was 
observed in the literature, while we find a weak crossover 
to a significantly higher value of ~ 0.43 eV for both 
translational and rotational correlation times. The for- 
mer value of the activation energy is consistent with re- 
sults from neutron scattering work^S whereas the latter 
is in agreement with activation energies determined in 
DS and NMR studieSf^^i^ The origin of these deviations 
is not yet understood. To clarify this point, it is desir- 
able to compare simulation results for various water and 
protein force fields in future work. 

In addition, we have demonstrated that an increase 
of the order of the HB network structure accompanies 
the change of the mechanism for the water dynamics and 
the weak crossover from fragile behavior to strong be- 
havior. Specifically, the percentage of water molecules 
exhibiting an ideal number of exactly 4 HB continuously 
increases upon cooling until it saturates at temperatures 
190 — 210 K, at which the temperature-dependent corre- 
lation times show the weak crossover. This agreement 
suggests an intimate relation between structure and dy- 
namics of the hydration water. When the temperature is 
decreased the energetic penalty a water molecule expe- 
riences when being in an imperfect HB environment be- 
comes more and more important, leading to constraints 
for the molecular positions and orientations. As a conse- 
quence, the mechanism for water dynamics changes from 
diffusive motion to jump motion and a local (3 relaxation 
sets in prior to the global a relaxation. These conclu- 
sions are in harmony with results of recent calculations 
for a cell model of waterjfiS. There, it was argued that, 
when the temperature is decreased, the number of HB 
increases and the orientational disorder decreases, result- 
ing in a rise of the activation energy and, hence, to non- 
Arrhenius behavior. When further cooling, the rate of 
change of orientational disorder reaches a maximum and, 
then, it rapidly drops to zero, leading to temperature in- 
dependent both disorder and activation energy at low 
temperatures. 

As compared to experimental results, the dynamics of 
SPC water is known to be too fast.^^'^'* Likewise, the wa- 
ter dynamics in our models is faster than that in hydrated 
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elastin and hydrated collagen. Specifically, t/j = 1 ns at 
200 K in the present study and at 250 K in our previous 

NMR spin-lattice relaxation analysis.— This discrep- 
ancy hampers straightforward comparison of computa- 
tional and experimental correlation times and crossover 
temperatures. Moreover, it is not clear whether results 
for the motional mechanism from MD simulations and 

NMR stimulated-echo experiments may be compared 
since water dynamics in the nanosecond regime and in 
the millisecond regime are studied, respectively. A com- 
mon activation energy of Ea « 0.45 eV suggests that both 
techniques probe the same low-temperature dynamical 
process of protein hydration water despite the different 
time windows. In harmony with the present results, we 
found in our previous NMR study that the mechanism 
for water dynamics changes from an isotropic motion to 
an anisotropic large-angle jump motion upon cooling.— 
While exact tt flips were ruled out as motional mecha- 



nism in the low-temperature Arrhenius regime, inexact 
TT flips may be consistent with our NMR data. More- 
over, exact TT flips cannot be observed in DS, in contrast 
to experimental flndings^ '^°i^^ but the small- angle wob- 
bling motion of the twofold symmetry axes of the wa- 
ter molecules associated with distorted tt flips would be 
probed by the latter method. Further experimental and 
computational studies are required to clarify these points. 
In particular, it is important to determine whether the 
existence of distorted tt flips is common to various estab- 
lished water and protein force fields. 
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